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Abstract: 


Nemipterus japonicus (Japanese threadfin bream) is one of the most significant fish species throughout the 
Indo-Pacific regions. No previous studies on the population dynamics of N. japonicus have been carried 
out in the Eastern Malaysia (Malaysian Borneo). The present study focused on the population dynamics 
of N. japonicus from April 2013 to March 2014 in the coastal area of Bintulu, South China Sea. Local 
fishermen provided samples, which were then analysed. The FiSAT II software was used to evaluate the 
collected length frequency data. The asymptotic length (L_) and growth coefficient (K) were 26.78 cm and 
0.85 yr-1, respectively. The growth performance index (’), total mortality co-efficient (Z), natural mortal- 
ity (M) and fishing mortality (F) were 2.785, 2.97, 1.63 and 1.34 per year, respectively. According to the 
estimates, the exploitation rate (E) for N. japonicus was 0.45, while the E,,, E, ,and E_, were determined 
as 0.72, 0.4 and 0.806, respectively. The obtained value of exploitation was less than E_, although slightly 
higher than E, ,. The recruitment characteristics of N. japonicus have been stable throughout the year, with 
one peak pulse. The results contribute to the knowledge of the dynamics of N. japonicus populations and 
are important for the sustainable management of the Sarawak fishing resources. 
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° from east Africa to the Indo-Malay Archipelago 
Introduction (ELHAWEET 2013, NING et al. 2015, IMTrIAz et al. 
The Japanese threadfin bream Nemipterus 2016, Fartvaretal. 2017, OGWANG etal. 2021). It 


japonicus (Bloch, 1791) (Nemipteridae) is char- 
acterised by pinkish body colour, silvery under- 
side, and 11-12 faint golden-yellow stripes run- 
ning down the body from below the head to the 
base (ELHAWEET 2013). Almost entirely of the 
Indo West Pacific waters is home to N. japonicus, 
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lives in shallow mud or sandy bottoms, at depths 
from 5 to 80 m. It is a commercially significant 
species in the South China Sea, where commer- 
cial trawling takes place in large numbers (Rus- 
SELL 1986, 1990, 1993, RUSSELL & Gouws 2013, 
NING et al. 2015). 
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Its importance has lead this species to be ex- 
tensively studied, although most of the stock and 
population assessment were carried out in the wa- 
ters of Indian subcontinent (MAsoomr et al. 2023; 
Raza et al. 2022; Rao et al. 2018; SEN et al. 2014; 
Jost 2010). Findings from research on population 
dynamics might help in planning and managing ma- 
rine resources by providing information on numer- 
ous population characteristics such as asymptotic 
length, growth coefficients, mortality rates from 
natural and fishery sources and the exploitation lev- 
el (NURUL AMIN et al. 1999, 2009, AL-MAMUN et 
al. 2021). More than 30 studies have been accom- 
panied on population dynamic of N. japonicus in 
different part of the globe, including the Northern 
Arabian Sea (IQBAL 1991), off Karnakata (ZACHARIA 
1998), off Visakhapatnam (RAJKUMAR et al. 2003), 
off Cochin (JosH1 2010), Madras (VIVEKANANDAN 
& JAMES 1986), Brunei (SILVESTRE & GARCES 2004); 
however, no studies from East Malaysia (Malaysian 
Borneo) have been carried out. The first report of 
this species in east Malaysian waters can be tracked 
back in 1977 (WEBER & JoTHY 1977). Later on, a 
little number of published documents dealt with dif- 
ferent aspects of N. japonicus in Sarawak or Sabah 
waters, including the status, breeding biology and 
genetic characteristics (KHATIB 2015, IMTIAz et al. 
2016, Lim et al. 2016). The breeding biology and 
feed preferences of N. japonicus were studied in 
Sarawak water (NETTELY et al. 2016, TONNIE et al. 
2018), while the genetic differentiation among the 
populations was investigated from Peninsular Ma- 
laysia (LIM et al. 2016). Population genetics research 
discovered similarities between N. japonicus found 
in the western Pacific and the Indian Oceans (IMTIAZ 
et al. 2016). However, the population structures of 
N. japonicus in east Malaysian waters, especially in 
Sarawak waters, are still unknown. 

The present study was performed in order to 
cover the gaps about the population dynamics of N. 
japonicus. One of the primary goals of this investiga- 
tion was to estimate the population characteristics and 
exploitation rate of N. japonicus with the purpose of 
evaluating the stock condition of N. japonicus in the 
coastal region of Bintulu, Sarawak, South China Sea. 


Materials and Methods 


Samples of threadfin breams N. japonicus were col- 
lected from commercial fish landing station at Bin- 
tulu (3°10’ 13.9” N, 113°02’22.9” E) and Kuala 
Nyalau village (3°38 26.9” N, 113°23’0.46” E) dur- 
ing the period of April 2013 to March 2014 (Fig. 1). 
The fishermen used a v-shaped purse seine net with 
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a mesh size of 4 inch. All samples were measured to 
the closest 0.1 cm for total length. Following that, 
one-centimetre intervals were utilised to partition 
the length frequency data into groups. 

As described by GAYANILO et al. (1996), the data 
were evaluated by means of the FISAT II in the com- 
puter software programme. In order to determine the 
asymptotic length (L_) and growth co-efficient (K) of 
the von Bertalanffy equation for increase in length, we 
used the ELEFAN-I software package included in the 
FiSAT II software package. The K scan process was 
carried out in order to determine whether an accurate 
assessment of the K value was obtained. The assessed 
L,, and K were accustomed to produce the ’ (growth 
performance index) with the formula, 


o’=2 log,,L_+ log, .K 


The possible longetivity (t,,.) of the investi- 
gated species was determined according to PAULY 
(1984) and TayLor (1958), who used the equation 
tax 3/K+ tO to arrive at his results. Finding the 
lengths at different ages was accomplished by the 
use of the inverse von Bertalanfy growth equation. 
The non-linear squares estimation approaches were 
accustomed to fit the VBGF in order to estimate the 
length-at-age curve. The following is the definition 
of the VBGF formula, 


L=L[l-e*“)] 


Where L, is the mean length at age t, while t is 
age of the species and t, is the hypothetical age at 
when the length is zero. 

The catch curve converted by the length was 
used to determine the total mortality coefficient (Z) 
(PauLy 1984). According to Pauly’s empirical re- 
lationship (PAULY & MARTOSUBROTO 1980), the M 
(natural mortality rate) was estimated with the em- 
pirical formula obtained from the empirical connec- 
tion as follows, 


Log,, M = -0.0066 — 0.279Log, ,L,, + 0.6543 Log,,K + 0.4634 Log, T 


In this equation, T is the annual mean tempera- 
ture (in °C) of the habitat water. 

As soon as Z and M values are known, it is 
possible to determine fishing mortality (F). F was 
calculated based on: 


F=Z-M 

The exploitation level (E) was estimated by the 
GULLAND (1971) relationship: 
E=F/Z=F/(F+M) 
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MALAYSIA 


Fig. 1. Sampling sites covered the coastal area of Bintulu in South China Sea. 


According to the approach of Paury & Sort- 
ANO (1986), the rising leftward arm of the length- 
converted catch curve was utilized to examine the 
likelihood of capture of each length class. By chart- 
ing the cumulative likelihood of capture against the 
mid-length at first capture (Lc), it was determined by 
corresponded to the cumulative probability at 50 %. 

A backward projection of cluster of acces- 
sible length frequency data on the length axis, as 
termed in software package, was used to identify 
the stock recruitment characteristics. The L,, K, and 
t0 (t=0) were the parameters used in the simulation. 
NORMSEP was used to determine the typical distri- 
bution of the recruitment configuration. 

The FiSAT protocol was used to perform the 
assessed length regulated virtual population analy- 
sis (VPA) and the cohort analysis for the assessed 
length regulated VPA. An analysis using the values 
of a (constant) and b (exponent) was performed on 
the values of L_, K, M, F, a (constant), and b (expo- 
nent). The value of t0 was assumed to be zero. 

According to the model developed by Beverton 
and Holt and updated by Pauy & SorRIANO (1986), the 
relative yield-per-recruit (Y/R) and relative biomass- 
per-recruit (B/R) were calculated and implemented 
into the FiSAT. It was necessary to provide the values 
of LC/L,, and M/K as input requirements for the op- 
eration to work. It was determined that the maximum 
permitted limit of exploitation (E,,,) was reached, 
resulting in the highest possible Y/R ratio. Addition- 


ally, the exploitation rates E,,, at which the marginal 


increase in relative Y/R is 10% of its value at E= 0 as 
well as E, ., at which the marginal increase in relative 
BYR is 50% of its value at E = 0, were calculated. 


Results 


Length frequency and population parameters 


The overall 1532 individual of N. japonicus were col- 
lected in present investigation. The lowest and highest 
total length were 10.5 cm and 25.5 cm respectively 
(Table 1). The asymptotic length (Loo) was 26.78, 
while the growth co-efficient (K) was 0.85 (Table 2). 


Growth parameters of Nemipterus japonicus 


According to the results, the observed and anticipated 
extreme lengths (L,,.) were 25.50 cm and 26.36 cm 
respectively (Fig. 2). The extreme length range was 
determined to be between 25.26 cm and 27.45 cm with 
a 95% confidence. The best estimate of the growth co- 
efficient (K) was 0.85 yr'. The response surface (Rn) 
was determined to be 0.240, and asymptotic length 
was found to be Lo=26.78 cm. During the current ex- 
periment, the computed value for growth performance 
index (’) of N. japonicus was 2.785 (Fig. 3). The re- 
constructed length-frequency histograms were super- 
imposed on top of the improved growth curve (Fig. 4). 


Mortality and exploitation rate 


The Z was found to be 2.97 yr'!, while the M was 
calculated to be 1.63 yr! (Fig. 5). The F was deter- 
mined to be 1.34 yr' based on the Z value. Based on 
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Fig. 2. The anticipated L,. value and the 95 % confidence level are determined by intersecting the inclusive L,. with 
the lines b and a, c, respectively. 


Table 1. The length frequency values of Nemipterus japonicus from the coastal area of Bintulu, Sarawak, South China Sea. 


Years 2013 2014 

Mid length (cm) Apr | May | Jun Jul Aug Sep Oct | Nov Dec Jan Feb Mar 
10.5 2 
11.5 11 3 
12.5 20 11 
13.5 3 17 1 3 3 7 
14.5 3 4 2 9 4 19 3 12 1 3 14 
15.5 6 5 6 22 2 18 24 10 15 7 15 
16.5 10 13 12 36 13 25 31 24 11 9 26 
17.5 13 18 18 29 19 24 21 29 21 2 25 30 
18.5 21 23 19 16 25 20 11 25 38 8 31 31 
19.5 25 24 19 4 18 18 8 12 19 6 25 34 
20.5 10 5 20 ) 20 11 8 4 11 16 15 17 
21.5 6 3 10 2 8 5 3 2 2 24 3 8 
22. 4 5 4 1 2 3 2 1 25 4 
23.5 3 2 1 1 2 1 12 2 
24.5 1 1 1 7 
25.5 1 1 1 

98 100 116 153 149 127 134 111 134 103 119 188 


Table 2. Various population characteristics of Nemipter- 1.00 
us japonicus collected from the coastal area of Bintulu, 0.90 | 
Sarawak, South China Sea. 0.80 
Parameters Values Ss sad i 
ie 26.78 % 0.60 | 
K 0.85 : 0:50 
Le 16.85 8 0.40 + 
o° 2.785 a 
Z 2.97 0.20 + 
M 1.63 ae | 
mae a 05 aE 5.0 10.0 
E 0.45 Growth constant K (/year) 
Em ws Fig. 3. The o’ of Nemipterus japonicus collected from 
N (ample size) 1532 Bintulu coast, Sarawak, South China Sea. 
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Fig. 4. The von Bertalanffy growth curves for Nemipterus japonicus (Loo = 26.78 cm and K = 0.85 yr') were overlaid 


on the reconstructed length frequency histograms. 


Length-Converted Catch Curve 
(for Z=2.97; M (at 27.0°C)=1.63; F=1.34; E=0.45) 
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Fig. 5. Least square linear regression was used to calcu- 
late the length converted catch curve; darker complete 
dots reflect those spots, whereas open dots show those not 
fully recruited or approaching Loo. Using our formulas, we 
came up with mortality rates M of 1.63 and Z of 2.97. 


these data, E of 0.45 was calculated, which was near 
to the optimal level (E=0.50). 


Length at first capture 


The L,, ,, =16.85 cm was derived from the examina- 
tion of the likelihood of capture, and this is the result 
that was acquired. In this study, the length at which 
75 % of N. japonicus remained in the gear was cal- 


culated to be L,,,, =17.86 cm (Fig. 6). 
Recruitment pattern 

The pattern of recruitment of N. japonicus was con- 
sistent throughout the year, with a single peak in Au- 
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Fig. 6. Probability of capture of each length class of the 
Nemipterus japonicus (L,, ,=15.68 cm, L,, =16.85 cm 
and L,, ,=17.86 cm) from Bintulu coast, Sarawak, South 
China Sea. 


gust? (Fig. 7). During the research period, the pro- 
portion of participants recruited ranged from 2.87 % 
to 21.32% of the total. The highest (21.32 %) was 
observed in November 2013 and the lowest (2.87 %) 
percent recruitment was in February 2014. 


Virtual population analysis (VPA) 


For the mid-lengths of 10.5 cm and 19.5 cm, the 
VPA of N. japonicus revealed that the lowest, F 
(minimum fishing mortality) and the F (maxi- 
mum fishing mortality) were reported at 0.004 and 
1.42 yr', respectively (Fig. 8). 
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Fig. 7. One peak pulse in the study year was seen in the 
characteristics of recruitment of Nemipterus japonicus in 
the Bintulu coast, Sarawak, South China Sea. 
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Fig. 8. Virtual population study of Nemipterus japonicus 
based on length in the Bintulu coast, Sarawak, South Chi- 
na Sea. 
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Fig. 9. The Y/R and B/R of Nemipterus japonicus by 
knife-edge procedure in the Bintulu coast, Sarawak, South 
China Sea. 
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Relative yield per recruit (Y/R) and biomass per 
recruit (B/R) 

The E_,, that provides the maximum relative yield 
per recruit was calculated to be 0.806, which relates 
to the maximum relative Y/R (Fig. 9). The maxi- 
mum and economic increases of Y/R at the exploita- 
tion rates were 0.806 (E_,.) and 0.72 (E, ,,), respec- 
tively. The exploitation rate that reduces the biomass 
to 50% of its unexploited (E, .) one was 0.4. 


Discussion 


During the study, 1532 specimens of N. japonicus 
were studied, with lengths ranging from 10.5—25.5 
cm. In accordance with the findings, the species of 
N. japonicus had L,, of 26.78 cm and K of 0.85 yr'. 
Previous research revealed that the values of L, var- 
ied from 24.5 cm to 34.5 cm and the values of K 
ranged from 0.27 year-1 to 1.0 year-1, respectively 
(Table 3). The lower value of L, was recorded in 
Bay of Bengal by Mustara (1994), while the higher 
was from Veraval water (SEN et al. 2014). For the 
K value the lower was from Pakistan water (KAL- 
HorO et al. 2014) and the highest was from Madras 
(VIVEKANANDAN & JAMES 1986). The value of L_ 
obtained from the calculation was 26.78 cm, and 
the value of K acquired from the computation was 
0.85 year!, which was considered the optimal value 
within the range for this species. 

The K termed ‘curvature parameter’, which 
regulates how rapidly the fish approaches its L, 
(RAHMAN & Cowx 2008) According to ADAMS 
(1980), species-specific ecological features, popula- 
tion size, and gene frequency all had a role in the 
disparities between documented L, and K values. 
These factors which influenced by habitat and natu- 
ral selection may cause to the different adaptation 
pattern during their life (ULASKI et al. 2020; CoNno- 
VER & SCHULTZ 1997). Recent studies suggested 
that, the L, and K of different fish species and bio- 
geographic regions also might influenced by the cli- 
mate change and global warming (VAN DENDEREN 
et al. 2020, WANG et al. 2020, BorserTTI et al. 2021). 

The recorded growth performance index of N. 
japonicus in the current study was high compared 
to that recorded in previous works (Table 3). Per- 
formance index is important to determine the fish to 
be considered for aquaculture suitability. Fish often 
lose weight but rarely lose length so, that the @ may 
be less liable to biases (MATHEWS & SAMUEL 1990). 
Because of the differences in size composition across 
species and the technique of study, the differences in 
growth parameter and $ within various locales may 
be linked to these differences (CARAMANTIN-SORI- 
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Table 3. The population parameters estimate of Nemipterus japonicus in various regions. 


SCS 


Reference Localities L,(em) | Kgr')| 6’ |M@r')| For’) | Zyr')| E 
KRISHNAMOORTHI (1971) Andhra-Orissa coast, India 20.90 0.65 | 2.45 - - - - 
KRISHNAMOORTHI (1971) Andhra-Orissa coast, India 30.30 0.25 | 2.43 - - - - 
KRISHNAMOORTHI (1973) India (Andra, Orissa) 30.50 0.31 | 2.74 - - - - 
LEE (1975) Hong Kong (females) 34.10 0.19 | 2.34 - - - - 
LEE (1975) Hong Kong (males) 38.20 0.13 | 2.28 - - - - 
Morty (1984) Kakinada 31.40 0.75 - - - - - 
PAULY & SANN AUNG (1984) Northern Myanmar (Burma) 37.00 0.24 | 2.51 - - - - 
PAULY & SANN AUNG (1984) Southern Myanmar (Burma) 37.00 0.24 | 2.52 - - - - 
INGLES & PAULY (1984) Manila Bay, Philippines 30.00 0.70 | 2.80 - - - - 
Corpuz et al. (1985) Carigara Bay, Philippines 23.50 | 0.730 | 2.61 - - - - 
Corpuz et al. (1985) Samar Sea, Philippines 26.50 0.60 | 2.62 - - - - 
Epwarbs et al. (1985) Gulf of Aden 29.10 0.31 | 2.42) 0.85 | 0.45 | 0.67 - 
VIVEKANANDAN & JAMES (1986) | Madras, India 30.50 1.00 | 2.97 | 2.50 0.45 | 2.98 - 
Morty (1987) Kakinada 33.90 0.52 | 2.78) 1.11 1.53. | 2.64 - 
DEVARAJ & GULATI (1988) India 29.80 0.82 | 2.86 - - - - 
Isa (1988) Kedah, West Malaysia 31.50 0.53 | 2.72) 1.18 | 2.84 | 4.02 | 0.71 
Isa (1988) Kedah, West Malaysia 31.40 0.55 | 2.73 | 1.21 2.51 | 3.72 | 0.68 
Humayun et al. (1989) Bay of Bengal, Bangladesh 26.50 0.60 | 2.62 | 1.32 3.39 | 4.71 - 
KHAN & MustaFa (1989) Bangladesh 24.16 1.06 | 2.79) 1.97 1.08 | 3.75 - 
IQBAL (1991) Northern Arabian Sea, Pakistan 28.80 0.46 | 2.58 - - - - 
GOPAL & VIVEKANANDAN (1991) | Veraval, India 33.70 0.73 | 2.91 - - - - 
MATHEWS & SAMUEL (1991) Kuwait 33.60 0.51 - - - 1.76 - 
Murty et al. (1992) Visakhapatnam, India 33.90 0.40 - - - - - 
BREIKAA (1992) Gulf of Suez 28.64 0.50 | 2.61 - - - - 
MustaFa (1994) Bay of Bengal 24.50 0.94 | 2.75) 0.78 | 0.55 1.33 | 0.41 
CHAKRABORTY (1995) Bombay, India 35.60 0.76 | 2.98) 1.55 | 2.03 | 3.58 - 
BREIKAA (1996) Gulf of Suez 29.27 0.46 | 2.60 - - - - 
LavaPIE-GONZALES et al. (1997) | Philippine 28.30 - - - - - - 
ZACHARIA (1998) es Seon amas 33.00 | 1.00 | - | 1.87 Z 5.65 | 0.68 
MustaFa (1999) Bay of Bengal, Bangladesh 25.60 0.94 | 2.79 - - - - 
AHMa~D et al. (2003) West Malaysia 34.80 0.85 - 1.11 2.41 | 3.52 - 
RasKumaR et al. (2003) Off Visakhapatnam 34.00 0.52 - 1.11 2.41 | 3.52 | 0.69 
EL-GANAINY & MEHANNA (2003) | Gulf of Suez 28.35 0.63 | 2.79 - - - - 
SILVESTRE & GARCES (2004) Brunei 28.50 0.65 | 2.72 | 1.37 0.31 1.68 | 0.18 
Josut (2010) Cochin, India (Male) 31.80 0.69 - 1.30 1.02 | 2.32 - 
Josut (2010) Cochin, India (Female) 26.50 0.77 - 1.30 | 0.76 | 2.06 - 
Josui (2010) Cochin, India (Pooled) - - - 1.30 1.87 3.35 - 
AMINE (2012) Gulf of Suez, Egypt 33.65 0.45 | 2.71) 0.53 1.22 | 1.75 | 0.69 
KaHoro et al. (2014) Pakistan waters 30.45 0.27 | 2.40 | 0.74 0.22 | 0.96 - 
SEN et al. (2014) Veraval, India 34.56 0.60 - 1.20 | 2.82 | 4.02 | 0.70 
Present study Pinmhncodsts Rest Maloyss 26.78 | 0.85 |2.79) 1.63 | 1.34 | 2.97 | 0.45 


Note: ‘-‘ (Not available) 
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ANO et al. 2008, MEHANNA et al. 2013, KALHORO et 
al. 2017, ZAIN & ABDULLAH 2019). 

For other species of Nemipterus such as N. de- 
lagoa, the estimated growth parameter for length 
frequency data at Tutocorin are L, =36.2 cm and 
K=1.06*". The growth parameters of N. mesop- 
rion from Mumbai waters (CHAKRABORTY 2002) 
based on data from 1989-1994 reported as L_=27.4 
mm, K=0.7628. The following parameter values of 
threadfin bream recorded from a Western Kalim- 
antan, South China Sea (PAULY & MARTOSUBROTO 
1980) was L, =24.5 cm and K = 0.42. When com- 
pared to the other related demersal species, the stock 
of N. marginatus evaluated can only be fished more 
tightly and the optimal fishing method, including 
the amount of effort expended, the mesh size and 
the types of gear employed (O’NEILL et al. 2018, 
Humpurigs et al. 2019, Lima et al. 2020, ROBERT et 
al. 2020, SALAZAR-PEREZ et al. 2020). Nemipterus 
randalli was studied along the Arabian Sea coast 
of Oman by AL-Kryumr et al. (2014) estimating that 
L,=22.12 cm and K=0.64 year'. The Z, as well as M 
and F, were assessed to be 3.0, 1.4, and 1.6 respec- 
tively. The exploitation ratio and rate were computed 
and were found to be 0.53 and 0.51, respectively, 
in this study. The L, of N. delagoa was the long- 
est measured among the various species, while the 
shortest measured length (L_) was N. randalli. 

The current research found that the F ranged 
from 0.22 to 2.82. The estimated M, which was com- 
parable with those reported in previous studies with 
the ranged between 0.74 and 2.52. The Z measured 
by length converted catch curve for N. japonicus in 
Bintulu coastal water was 2.97 yr! and comparable 
to the value at Madras (2.985 yr!) reported by VIVE- 
KANANDAN & JAMES (1986). The exploitation level 
of 0.45 shows that the V. japonicus fishery in Bintu- 
lu coastal waters was being under exploitation dur- 
ing the study periods. This was based on the premise 
that the stock was maximally exploited when the F 
was equal to M, or when E= (F/Z) = 0.5, or when F 
was equal to M. 

Variation in M may be interpreted as a natural 
phenomenon influenced by both density-dependent 
variables such as predation and food availability as 
well as density-independent causes such as diseases 
and natural catastrophes, among others. It was ob- 
served that the mortalities varied within same spe- 
cies in different location (SEN et al. 2014). 

A single recruitment event each year has been 
seen for N. japonicus (Fig. 6), which is in keeping 
with previous research by SEN et al. (2014) from Ve- 
raval water which found that the recruitment pattern 
for this species occurs constantly year round, with a 
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top peak in May to July. It was found that between 
2.71 and 16.25 % of those surveyed were success- 
ful in getting jobs. The month of June had the great- 
est rate of recruiting (16.25 %), while the months of 
January and November had the lowest rates. 

Among other things, length at first capture is 
a significant input in the evaluation of relative Y/R 
and comparable B/R when using the knife edge se- 
lection procedure because it is assumed that fishes 
of smaller size than length at first capture were not 
caught by fishing gears (BEDDINGTON & KIRKWOOD 
2005, CARRUTHERS et al. 2014, UpoH & UKpaTu 
2017). 

It was determined in the present study that the 
greatest permitted limit of E, that would result 
in the highest possible relative Y/R was assessed 
to be 0.8006. This is in contrast to the E that was 
found to be 0.45. Using the ideal value of 0.50 as a 
benchmark, it can be determined that the fisheries 
was being exploited optimally. This conclusion was 
reached after conducting a mortality-based exploi- 
tation rate study and calculating the relative Y/R. 
The results of these analyses suggested that the this 
species are exploiting below its ideal level. As a re- 
sult, it is recommended that the fishing pressure on 
the stock be maintained at the present level or re- 
duced. Furthermore, fish were captured before they 
reached maturity size, which may have an impact on 
the future population and stock size (HUTCHINGS & 
REYNOLDS 2004, PINsky et al. 2011, BEN-HASAN et 
al. 2021). Because of the decline in the mesh size of 
trawl nets at the cod-end, it is possible to see a fall in 
the lengths at first catch (LC,,), which may lead to 
recruitment overfishing over an extended period of 
time (O’NEILL & KYNOCH 1996, HERRMANN 2005, 
SALA et al. 2016, ROBERT et al. 2020). 


Conclusions 


Findings from the present study indicate that N. 
japonicus particularly peak during the post-mon- 
soon period. The exploitation rate for the species 
was 0.45, indicating that the resource is being under 
exploitation. It is feasible to get further information 
on the other Nemipterus spp. of Sarawak using SCS 
in order to ensure their sustainable harvest and to 
estimate the probable stock. 
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